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Abstract 

The last decade has seen the emergence of a wide 
range of automated calculations for supersymmetric 
extensions of the Standard Model. This guide con- 
tains a brief summary of these, with the main focus 
on hadron collider phenomenology, as well as a brief 
introduction to the so-called SUSY Les Houches Ac- 
cord. See also the Les Houches Web Repository for 
BSM Tools: 

http : //www. ippp . dur . ac .uk/montecarlo/BSM/ 

1 Introduction 

Among the most enticing possibilities for observable 
New Physics both at the Tevatron and at the Large 
Hadron Collider is supersymmetry (SUSY); for re- 
views, see e.g. [1-3]. At the most fundamental level, 
imposing supersymmetry on a quantum field theory 
represents the most general (and only) possible way 
of extending the Poincare group of space-time sym- 
metries [4,5], at the same time as it furnishes a de- 
sirable relation between the bosonic and fermionic 
degrees of freedom. Empirically, however, supersym- 
metry can at most be a broken symmetry if it ex- 
ists in Nature, due to the non-observation of a mass- 
degenerate (or lighter) spin-partner for each of the 
Standard Model (SM) particles. 

However, even a softly broken supersymmetry can 
have quite amazing properties, as long as the mass 
splittings introduced by the breaking are smaller than 
a TeV or so. Among the most well-known conse- 
quences of such a type of supersymmetry are radia- 
tive breaking of electroweak symmetry, an elegant so- 
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lution to the so-called hierarchy problem, a natural 
weakly interacting dark matter candidate (in theo- 
ries with conserved i?-parity), and unification of the 
strong, weak, and electromagnetic gauge couplings at 
a (very) high energy scale. 



For collider phenomenology, the most immediately 
relevant consequences are 1) an extension of the Stan- 
dard Model Higgs sector to (at least) 2 doublets, 2) 
promotion of each of the Standard Model fields (plus 
the extra Higgs content) to superfields, resulting in a 
spin-partner for each SM particle, with mass split- 
tings inside each boson-fermion doublet 1 TeV, 
and 3) the special properties which accompany a con- 
served i?-parity, namely production of the new states 
only by the pair, followed by individual cascade de- 
cays down to the Lightest Supersymmetric Particle 
(LSP) which is stable and (usually) escapes detec- 
tion. 



The large interest in (N = 1) supersymmetric ex- 
tensions of the SM and their phenomenological con- 
sequences has carried with it the need for automated 
tools to calculate supersymmetric mass and coupling 
spectra, cross sections, decay rates, dark matter relic 
densities, event rates, precision observables, etc. To 
handle the cross-communication between the many 
tools, the so-called SUSY Les Houches Accord [6-8] 
(SLHA) is now in widespread use. Section contains 
a brief introduction to this Accord. Next, in Section 
13 an overview of the presently available state-of-the- 
art tools is given, divided into four main categories. 
A more extensive collection of tools for BSM physics 
as well as an online repository can be found in [9]. 
Another recent and comprehensive tools review is the 
Les Houches Guidebook to MC Generators [10]. 
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2 The SUSY Les Houches Ac- 
cord 

Given the long history of the subject, it is not sur- 
prising that several different conventions for defining 
supersymmetric theories have been proposed over the 
years, many of which are in active use by different 
groups. While this is not a problem per se (unique 
translations can always be constructed), it does en- 
tail a few practical problems, in particular when the 
results of one group are compared to (or used for) the 
calculations of a different group. 

In addition, even when the theoretical conventions 
are identical, there remains the purely technical is- 
sue that each program has its own native way of in- 
putting and outputting parameters, each of which is 
unintelligible to most other programs. 

The SLHA was proposed to solve both these prob- 
lems. Due to the large parameter space of uncon- 
strained supersymmetric models, the SLHA in its 
present form [6] (SLHA1) is limited to the MSSM, 
with conservation of i?-parity, CP, and flavour. Ex- 
tensions to more general models are underway [8] 
(SLHA2). 

Technically, the Accord is structured into 3 ASCII 
files (or strings): 1) model definition and measured 
SM parameters, 2) SUSY mass and coupling spec- 
trum, and 3) decay tables. Though admittedly not 
elegant, the ASCII format was chosen for its robust- 
ness across platforms and compilers. In general, all 
input parameters used for a calculation are copied to 
the output, so that any subsequent calculation also 
has access to the exact input parameters used for the 
previous one. 

2.1 The SLHA Conventions 

The backbone of the Accord is a unique set of con- 
ventions for defining the supersymmetric parameters, 
fields, and couplings. These conventions, which have 
also been adapted for the so-called SPA project [11], 
largely resemble the widely used Gunion-Haber con- 
ventions [12], with a few differences as noted explic- 
itly in [6]. Simply stated, to define a SUSY model, 
one needs the field content, the Superpotential, the 
SUSY breaking terms, and the gauge couplings. For 



the field content, the SLHA assumes that of the 
MSSM, while SLHA2 will include extensions for the 
NMSSM. 

The MSSM Superpotential is specified by the mea- 
sured SM particle masses (giving the Yukawa cou- 
plings) and by the fi term. At present, only the third 
generation Yukawas arc included. The gauge cou- 
plings are specified in terms of Mz, Gf, a s (Mz) MS , 
and the fine structure constant at zero momentum 
transfer. All of these are the standard SM ones that 
one can get from a review text, i.e. no SUSY correc- 
tions should be included here. SLHA2 will include 
masses for all 3 generations, as well as the CKM ma- 
trix. 

The SUSY breaking terms can either be specified 
by giving the parameters for a minimal version of a 
particular SUSY breaking model (SUGRA, GMSB, 
or AMSB), or individually, either by starting from a 
minimal model and successively adding non-universal 
terms or by simply giving all terms explicitly. For 
higher-order calculations, these parameters are inter- 
preted as given in the modified dimensional reduction 
(DR) scheme [13-15], either at the (derived) unifica- 
tion scale or at a user-specifiable scale. As mentioned, 
CP, i?-parity, and flavour are assumed conserved in 
SLHA1. 

In the spectrum output, three kinds of parameters 
are given: 1) pole masses of all (s)particles, 2) mixing 
matrices, and 3) Lagrangian parameters. While the 
precise definition of the mixing matrix elements are 
left up to each spectrum calculator, the Lagrangian 
parameters are defined as DR ones at one or several 
user-specifiable scales Q. 

2.2 The SLHA Decay Tables 

A somewhat separate and self-contained aspect of the 
SLHA is the possibility to pass total widths and par- 
tial branching ratios via a file structure similar to 
that of the rest of the Accord. A common use for 
this is to improve or extend the width calculations of 
an event generator by the numbers calculated with a 
specialised package. 

Note! An important potential pitfall when using 
these files is on-shell intermediate resonances in fi- 
nal states with more than 2 particles. If not treated 
properly, large problems both with double-counting 
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and with incorrect population of phase space can oc- 
cur. Please see [6] for an explicit description of the 
correct procedure to adopt in these cases. 

3 Computing SUSY 

This Section contains an overview of SUSY calcula- 
tional tools, divided into 1) spectrum calculators, 2) 
observables calculators, 3) matrix element and event 
generators, and 4) fitting programs. For links and 
references, the reader should consult the recently con- 
structed online repository for BSM tools [9]. 

3.1 Spectra 

Given assumptions about the underlying supersym- 
metric theory (field content, superpotential, super- 
symmetry breaking terms) and a set of measured pa- 
rameters (SM particle masses, gauge couplings, elec- 
troweak symmetry breaking), the masses and cou- 
plings of all particles in the spectrum can be com- 
puted. This is the task of spectrum calculators, also 
called RGE packages. 

The most commonly used all-purpose spectrum 
calculators are Isajet [16], SoftSusy [17], SPheno 
[18], and SuSpect [19], all compatible with SLHA. 
In general, the codes agree with each other to within 
a percent or so, though larger discrepancies can occur 
in particular at large tan /?. For mSUGRA, a useful 
online tool for comparison between them (and dif- 
ferent versions of them) exists [20] . Recent detailed 
comparison studies can be found in [21-23]. Though 
Pythia also contains an internal spectrum calcula- 
tor [24], the resulting spectrum is very approximate 
and should not be used for serious studies. 

There are also a few spectrum calculators with 
more specialised areas of application, here CPSu- 
perH [25], FeynHiggs [26], and NMHDecay [27]. 
NMHDecay computes the entire mass spectrum in 
the NMSSM (and has a limit which is equivalent to 
the MSSM), but couplings and decay widths are so 
far only calculated for the Higgs sector, though im- 
provements are underway. NMHDecay is compat- 
ible with an extension of the SLHA [8]. The pro- 
gram FeynHiggs deals with the Higgs sector of the 
MSSM, for which it contains higher precision calcula- 



tions than the general-purpose programs mentioned 
above. It is also able to handle both minimal flavour 
violation (MFV) and CP violation, and is compatible 
with the SLHA, hence can e.g. be used to provide a 
final adjustment to the Higgs sector of a general spec- 
trum calculated by one of the other codes. Finally, 
CPSuperH deals with the Higgs sector in the MSSM 
with explicit CP violation and contains a number of 
refinements which makes it interesting also in the CP 
conserving case. 

3.2 Observables 

Programs that calculate one or more of the follow- 
ing: cross sections, decay partial widths, dark mat- 
ter relic density, and indirect /precision observables. 
Note that we here focus on calculations relevant for 
hadron colliders and that matrix element and event 
generators, which also calculate many of these things, 
are treated separately below. 

For hadron collider cross sections, Prospino [28] 
can be used to calculate inclusive SUSY-NLO cross 
sections, both total and differential. It also calculates 
the LO cross section and gives the corresponding K- 
factor. 

For decay partial widths, several specialised pack- 
ages exist. For the MSSM, SPheno calculates tree- 
level decays of all (s)particles (soon to include RPV 1 ), 
S Decay [29] computes sparticle decay widths in- 
cluding NLO SUSY-QCD effects, and both Feyn- 
Higgs [26] and HDecay [30] compute Higgs par- 
tial widths with higher-order corrections. NMHDe- 
cay [27] computes partial widths for all Higgs bosons 
in the NMSSM. 

For the density of dark matter, DarkSUSY [31], 
IsaTools [32], and MicrOMEGAs [33] represent 
the publically available statc-of-thc-art tools. All of 
these work for the MSSM, though a special effort has 
been put into MicrOMEGAs to make it easily ex- 
tendable [34] , recently resulting in an implementation 
of the NMSSM [35], and work on CP violation is in 
progress. 

For precision observables, MicrOMEGAs includes 
calculations of (g — 2)^, b — ► S7, B s — > and 

1 RPV in SPheno is not yet public, but a private version is 
available from the author 
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cross-sections for ncutralino annihilation at v^O, rel- 
evant for indirect detection of neutralinos. NMHDe- 
CAY includes a check against LEP Higgs searches, 
b — > S7, and can be interfaced to MicrOMEGAs 
for the relic density. Isajet/IsaTools include cal- 
culations of b — > S7, (g — 2) p , B s — > — ► 
t + t _ , and neutralino-nucleon scattering cross sec- 
tions. SPheno includes b — > S7, (g — 2) M , as well as 
the SUSY contributions to the p parameter due to 
sfcrmions. FeynHiggs also evaluates the contribu- 
tion to electroweak precision observables via Ap as 
well as (g — 2)„ (with two- loop corrections). Finally, 
SuSpect also includes a calculation of b — > S7. 

3.3 Matrix Element and Event Gen- 
erators 

By a matrix element generator, we here understand a 
program that, given a set of fields and a Lagrangian, 
is able to generate Feynman diagrams for any pro- 
cess and square them. Note, however, that many of 
the codes are able to do quite a bit more than that. 
An event generator is a program that, given a matrix 
element, is able to generate a series of random exclu- 
sive events in phase space, often including resonance 
decays, parton showers, underlying event, hadronisa- 
tion, and hadron decays. 

The automated tools for generating matrix ele- 
ments for SUSY are Amegic++ [36], CalcHEP 
[37], CompHEP [38], Grace-SUSY [39], SUSY- 
MadGraph [40], and O'Mega [41]. All of these 
work at Leading Order, except Grace, and all cur- 
rently only deal with the MSSM, except CalcHEP 
which contains an NMSSM implementation. 

CalcHEP and CompHEP provide internal event 
generators, while the event generator Sherpa [42] 
is built on Amegic++, GrOppa [43] builds on 
Grace, MadEvent [44] builds on MadGraph 
(work is in progress to extend this to SUSY- 
MadGraph), and Whizard [45] builds on O'Mega. 
Of these, most are matrix-element-level event gener- 
ators, that is they provide events consisting of just a 
few partons and their four-momenta, corresponding 
to the given matrix element convoluted with phase 
space. These events must then be externally in- 
terfaced [46] e.g. to Pythia or Herwig for reso- 



nance decays, parton showering, underlying event, 
and hadronisation. The exception is Sherpa, which 
contains its own parton showers and underlying-event 
models (similar to the Pythia ones), and for which 
a cluster-based hadronisation model is being devel- 
oped. 

In addition, both Herwig [47] and Pythia [48] 
contain a large number of internal hardcoded leading 
order matrix elements, including i?-parity violating 
(RP V) decays in both cases [49-53] , and RPV single 
sparticle production in Herwig [49]. In Pythia, the 
parton shower off SUSY resonance decays is merged 
to the real NLO jet emission matrix elements [54], 
an interface to CalcHEP and NMHDecay exists 
for the NMSSM [55], and an implementation of the 
hadronisation of i?-hadrons is available [56,57]. 

Two other event generators should be mentioned. 
Isajet [16] also contains a large amount of SUSY 
phenomenology, but its parton shower and hadroni- 
sation machinery are much less sophisticated than 
those of Herwig, Pythia, and Sherpa. The active 
development of Susygen [58] (which among other 
things includes RPV single sparticle production) is 
currently at a standstill, though basic maintenance 
is still being carried out. 

3.4 Fitters 

Roughly speaking, the tools described above all have 
one thing in common: given a set of fundamental 
parameters (themselves not directly observable) they 
calculate the (observable) phenomenological and ex- 
perimental consequences. However, if SUSY is at 
some point discovered, a somewhat complementary 
game will ensue: given a set of observed masses, cross 
sections, and branching ratios, how much can we say 
about the fundamental parameters? 

The fitting programs Fitting- [59] and Sfitter 
[60] attempt to address this question. In a spirit 
similar to codes like Zfitter [61], they combine the 
above tools in an automated statistical analysis, tak- 
ing as input a set of measured observables and yield- 
ing as output a set of fundamental parameters. 

Obviously, the main difficulty does not lie in deter- 
mining the actual central values of the parameters, 
although this can require significant computing re- 
sources in itself; by far the most important aspect 
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of these tools is the error analysis. Statistical uncer- 
tainties can be treated rigorously, and are included 
in both programs. Theoretical and systematic uncer- 
tainties are more tricky. In a conventional analysis, 
these uncertainties are evaluated by careful consid- 
eration of both the experimental setup, and of the 
particular theoretical calculations involved. In an au- 
tomated analysis, which has to deal simultaneously 
with the entire parameter space of supersymmctry, 
a 'correct' evaluation of these errors poses a truly 
formidable challenge, one that cannot be considered 
fully dealt with yet. 
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